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ABSTRACT 


When  the  cost  of  reducing  the  duration  of  activities  is  convex 
and  nonlinear,  it  may  be  advisable  (to  reduce  the  computing  burden) 
to  seek  a  "sat Isfising"  answer,  in  which  the  project  is  "compressed" 
to  a  desired  completion  time  with  prespecified  tolerable  relative 
error.  We  treat  the  problem  of  constructing  the  optimal  first  degree 
interpolating  linear  spline  that  guarantees  such  maximal  error,  and 
consider  various  possible  refinements. 
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INTRODUCTION 


A  significant  recent  trend  in  the  study  of  activity  networks 
(ANs)  is  the  shifting  of  focus  from  insisting  on  achieving  the  opti¬ 
mum  to  requiring  only  approximations  that  deviate  from  the  exact 
results  by  a  known  (perhaps  predetermined)  error;  i.e.,  satisfising 
instead  of  optimizing.  In  the  overwhelming  majority  of  cases, 
asking  for  an  approximation  is  more  realistic  than  demanding  the 
exact  answer  because  of  the  very  approximate  nature  of  the  input 
data. 

It  is  a  truism  that  project  compression  under  linear  or  piece- 
wise  linear  cost  functions  of  the  form  C(y)  =  b  -  ay;  £  <  y  <  u; 
a,b  >  0  is  considerably  easier  to  resolve  optimally  than  under 
nonlinear  cost  functions.  The  analysis  of  [  2 ]  and  [3],  carried 
under  the  simplifying  assumption  of  quadratic  cost  functions,  should 
amply  demonstrate  this  fact,  if  such  demonstration  were  needed!  The 
natural  question  then  is:  What  if  C(y)  is  not  quadratic,  though 
still  convex  decreasing  as  y  increases  from  i  to  u?  For  instance, 
suppose 
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C(y)  =  a/ (b  +  ry)  ;  0  <  H  <  y  <  u  <  «  ;  a,b,r  >  n  (1) 

What  can  be  said  about  the  optimum  in  this  case? 

One  may  wish  to  persist  in  applying  the  theoretical  con¬ 
structs  of  [ 2 ]  and  [3],  which  are  indeed  applicable  in  toto. 
Unfortunately,  the  Kilter  Diagram  (which  is  simply  the  plot  of  the 
derivative,  -dC/dy  versus  the  duration  y)  will  now  possess  a  nonlinear 
segment  in  region  R,  as  shown  in  Fig.  1  for  C(y)  of  (1),  which  would 
necessitate  the  solution  of  a  system  of  nonlinear  equations  in  the 


Figure  1  approximately  here 

"flows"  {f^ j }  at  each  iteration;  an  onerous  task  at  best. 

The  other  alternative  is  to  approximate  the  cost  function 
C(y)  by  a  continuous  piecewise  linear  and  convex  function  (i.e., 
linear  spline)  that  is  "optimal"  in  some  sense.  There  are  two 
immediate  questions  that  present  themselves:  The  first  is  to 
define  the  criterion  of  optimality  of  the  approximation,  and  the 
second  is  to  define  the  sense  of  the  approximation  itself. 

This  paper  answers  these  two  questions,  then  proceeds  to  demonstrate 
the  application  of  these  answers  to  achieve  the  desired  approxima¬ 
tion. 

Although  our  discussion  is  couched  throughout  in  the  ver¬ 
nacular  of  ANs,  the  approach  has  broad  applicability  to  a  wide  class 
of  convex  separable  programming  problems. 


A 
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THE  CRITERION 

Let  C(y)  be  a  nonlinear  convex  decreasing  function  on  [*.,  u]. 
Let  the  points  y^y^  •  •  .  »yn+1  satisfy  i  =  yQ  <  y1  <  .  .  .  <  yn 
<  yn+1  =  u  and  let  c^,  1  <  i  <  n  +  1  be  arbitrary  real  numbers. 
Define 

h(y)  =  cilL  +  (ci  -  ci_1)(y  -  yi_1)/(yi  -  y^)  (2) 

for  ye[yil,yi]  ,  1  i  i  i  n  +  1 

Then  h  is  a  continuous  piecewise  linear  function  called  a  first 

degree  spline  with  knots  {y^},  0  <_  i  <_  n  +  1  that  joins  the  points 

(y. ,  c.),  0  <  i  <  n  +  1.  Since  the  knots  yn  and  y  ,  are  fixed,  the 

problem  is  to  choose  the  (variable)  knots  yyY2>  •  •  •  »yn  31,(1  the 

ordinates  •  •  •  >cn+i  so  tliat  the  corresponding  spline  h  is 

as  close  to  C  as  possible  in  the  sense  of  the  uniform  norm;  that  is, 

we  choose  h  to  minimize  the  maximum  error: 

min:  e  =  max  | C (y)  -  h(y)| 

yz[l,  u] 

=  II  C  -  h  ||  (3) 

where  h  is  a  first  degree  spline  with  n-1  variable  knots  in  the 
interval  [l,  u].  If  such  an  h  exists,  it  is  called  a  best 
Chebychev  (uniform,  minimax)  variable  knot  first  degree  spline  (FDS) 
approximation  to  C  on  [l,  u].  The  value  of  t  is  called  the  error 
or  deviation  of  the  approximation  h. 
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Traditionally  in  spline  approximation  the  values  at  the  knots 

satisfy 

=  C(y^)  ,  0  <  i  <  n  +  1  (4) 

In  this  case,  h  is  called  a  first  degree  interpolating  spline  (FDIS). 
If  a  spline  li  exists  which  minimizes  (3)  and  satisfies  (4),  we  call 
h  a  best  Chebychev  variable  knot  interpolating  spline  approximation 
to  C  on  [l,  u] .  If  we  relax  the  restriction  that  c^  =  C(y^), 

0  <  i  <_  n  +  1,  then  we  call  h  a  best  Chebychev  first  degree  approxi¬ 
mating  spline  (FDAS)  with  variable  knots.  Three  important  facts; 
that  are  known  about  optimal  Chebychev  splines  are  worth  recalling 
because  of  their  relevance  to  our  subsequent  development: 

1.  The  errors  in  all  segments  of  the  approximation  are 
equal  [  5 ] . 

2.  The  maximum  error  in  FDIS  is  exactly  twice  that  in  FDAS, 
which  yields  the  least  error. 

3.  For  purposes  of  optimization,  we  may  use  any  linear  spline 
approximation  that  differs  from  the  FDAS  by  a  constant 
(the  so-called  "simple  linear  translation")  [4]. 

Consequently  we  shall  deal  exclusively  with  FDIS.  This  choice 
will  become  even  the  more  advisable  as  our  discussion  progresses. 

In  ref.  [6]  Salem  gives  elementary  constructions,  both 
graphical  and  analytical,  for  the  determination  of  the  optimal  FDIS 
(and  FDAS)  that  guarantee  an  error  of  value  e.  (As  a  consequence, 
the  number  of  knots  is  variable.)  For  our  purposes,  we  briefly  state 
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the  concepts  underlying  his  graphical  approach  to  the  determination 
of  the  FDIS.  The  reader  should  have  no  difficulty  in  constructing  the 
analytical  procedure,  which  follows  the  graphical  construction  rather 
faithfully,  and  is  stated  in  full  detail  in  [6].  (In  fact,  we  do 
just  that  for  the  specialized  cost  function  in  the  example  below.) 

Suppose  that  C(y)  is  plotted  to  some  reasonable  scale,  as 
shown  in  Fig.  2.  From  the  left-most  point  (£;  C (£)),  drop  an  e  and 


Figure  2  approximately  here 

draw  a  tangent  to  the  function  C(y);  and  /rum  the  point  (A;  C(J.)), 
draw  a  parallel  line  to  that  tangent  to  intersect  C(y)  at  the  point 
(y1,C(y1)) .  This  is  the  first  segment  of  the  FDIS.  The  construc¬ 
tion  is  now  repeated  starting  with  the  point  (y^CCy1));  and  is 
stopped  when  the  upper  limit  u  is  reached.  This  "last  segment"  of 
the  linear  spline  bears  close  scrutiny,  since  it  may  lead  to  further 
improvement  in  the  approximation.  We  shall  do  this  later,  but  for 
the  moment  it  should  be  evident  that  the  above  construction  yields 
a  FDIS  whose  maximum  deviation  from  C(y)  is  at  most  e.  (The  maximum 
deviation  in  the  last  segment  may  be  strictly  <  r.) 

THE  SENSE  OF  APPROXIMATION 

Our  concern  is  next  focused  on  defining  the  sense  of  approxi¬ 
mation  that  is  appropriate  for  our  application,  and  on  determining 
the  values  e^j;(ij)eA,  that  realize  it. 

Recall  that  the  problem  of  optimal  project  compression  may 


be  stated  mathematically  as  (see,  e.g.,  ref.  [l]): 
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Program  P:  Minimize  C(y)  = 

(ij)eA 

C.  .  (y.  .) 

13  U13 

(r>) 

Sub j  ect  to  t .  -  t  •  + 

i  3 

yii  = 

0  ;  (ij)eA 

(6.1) 

-  t,  +  t 

1  n 

= 

T 

s 

(6.2) 

H  .  .  <  V .  . 
I]  - 

<  u.  . 

-  H 

(6.3) 

Because  of  the  difficulty  in  a  frontal  attack  on  Program  P  because 
of  the  nonlinearity  of  the  objective  function,  we  propose  to  approach 
it  through  approximating  each  individual  C^.(y.j)  by  a  FDIS,  denoted 
by  h^Cy  —  ).  Thus  we  define  the  following  surrogate  program: 

V 

Program  S:  Minimize  H(Y)  =  ^  h..(y..)  (7) 

(ij)eA  1:1  ” 

Subject  to:  (6.1)  -  (6.3) 

Let  v(P)  and  v(S)  denote  the  minima  of  their  respective  programs. 
There  are  two  possible  criteria  that  measure  the  "goodness"  of  this 
approximation: 


Criterion  1: 


Find  (ij)eA,  such  that  the  absolute  deviation  from 

the  exact  optimum  is  bounded  by  6  : 


|v(S)  -  v(P)  )  jc  6  ;  <5  >  0 


(8) 


Criterion  2: 


Find  e. (ij)eA,  such  that  the  relative  deviation  from 
the  exact  optimum  is  bounded  by w  : 


|v(S)  -  v(P)|/|v(P)|  <  «  ;  w>0  ,  |v(P) |  >  0 


O') 
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In  (8),  the  scalar  6  is  an  "absolute"  measure  since  it  lias 
the  same  dimension  as  the  function  C  (say  dollars),  while  the  scalar 
<0  of  (9)  is  dimensionless  since  it  represents  the  ratio  of  the  (abso¬ 
lute)  error  of  approximation  to  the  minimal  value  of  program  P; 
typically  we [.01,  .10],  Clearly,  Criterion  1  is  useful  when  the 
analyst  has  a  fairly  precise  idea  about  the  value  of  v(P)  on  which 
he  can  base  the  value  of  d.  Criterion  2,  on  the  other  hand,  is  free 
of  such  prior  knowledge,  since  it  requires  only  the  specification  of 
a  tolerable  relative  deviation  from  the  optimum  of  program  P  (a 
" satisfising"  relative  error)  whatever  that  value  may  be. 

Consider  Criterion  1:  We  wish  to  find  e. -  such  that 

il 

|v(S)  -  v(P) |  =  v(S)  -  v(P)  <  6 


The  reader  will  now  realize  the  added  rationale  for  choosing  the 

FDIS,  which  enables  us  to  replace  the  absolute  difference  by  the 

straight  difference,  since  h. . (y. .)  >  C. . (y. .)  for  £. .  <  y. .  <  u. . . 
h  ij  ij  -  il  —11  —  il 


v(S)  = 

Y  and  Y* 


Z_.  h .  -  (y .  . )  >  X!c..(y..)  >  X 
(ij)  13  13  "  it  13  13  "  Cij) 


C. . (y. .)  =  v(P),  where 
ij'ij 


are  the  optimizing  vectors  for  Programs  S  and  P,  respectively. 


But 


v(S)  -  v(P) 


min  ^C.-Cy..) 

yen  (ij)  yen  (ij) 


where  n  is  the  space  of  feasible  solutions  defined  by  (6.1)  -  (6.3). 
If  the  optimum  of  program  S  is  realized  at  ?,  and  the  optimum  of 


8 


* 

program  P  is  realized  at  Y  ,  then  for  each  activity  (ij)  the  contri- 

*  * 

bution  to  the  difference  is  h.  . (y..)  -  C..(y. which  is  obviously 

13  i3  ij  7iJ 

<_e^.  Consequently,  we  conclude  that 

|v(S)  -  v(P)  |  i  I  Xi  eijl  -  X  ^ij^ 

(ij)  (ij) 

Therefore,  by  taking  =  <5/m  inequality  (8)  will  be  satisfied, 
where  m  is  the  number  of  activities  in  the  network,  m  =  [ A [ . 

Next,  consider  Criterion  2,  which  is  the  one  adopted  hence¬ 
forth.  3y  construction  of  the  FDIS, 


h.  .  (y.  .)  -  C.  .  (y .  .)  <  e.  . 
i)  l]  13  13  “  13 


Summing  overall  (ij)eA, 


V1  V1  V 

/  i  h..(y.O  <  /,C..(y..)  +  /__i  e... 

'  13  13  -  1  13  13  13 


E  .  . 

13 


Minimizing  both  sides  of  the  inequality  results  in 
v(S)  _<v(P)  +  X! 
which  we  rewrite  as 

y 

<  1  +  <_i  e .  .  / 


(10) 


v (S) /v (P) 


13 


/v(P)  ;  v(P)  >  0 


(10') 


Actually,  should  be  written  (e^j|u),  and  v(S)  should  be  written 

(v(S)  |  {e^ }) ,  since  is  indeed  a  function  of  to,  and  v(S)  is  a 
function  of  the  errors  {e,..}.  We  forfeit  such  rigor  of  expression 
for  the  sake  of  clarity  of  exposition. 
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From  the  statement  of  Criterion  2,  we  have 


v(S)/v(P)  <_  1  +  u  (21) 

Equating  the  r.s.  of  (10)  and  (11),  and  assuming  for  simplicity  (see 
below)  that  e^  =  e,  the  same  value  for  all  activities,  we  conclude 
that 


e  <  oj.v(P)/ra  ;  m  =  j A J 


Since  the  value  of  v(P)  is  not  known,  any  lower  bound  on  that  value 
will  yield  a  more  "conservative"  e.  Such  a  lower  bound  is  immediately 


available  in 


/  C. . (u. . ) . 
Z_i  ij  n  ’ 

(ij) 


and  we  put 


c  =  u.C(U)/ m  <  w.v(P)/m 


(12) 


where  C(U)  =  C^(in_.).  This  is  the  e  to  be  used  in  the  deter - 

(ij) 

mination  of  the  Chebychev-optimal  FDIS. 

A 

Once  the  approximate  problem  S  is  solved  and  the  values  Y 
and  v(S)  are  in  hand,  one  may  obtain  the  a  posteriori  bounds 


v(S)  - 


Z  eij  ^v(p)  ±c 


C(Y  ) 


(IS) 


(ij) 


The  right  inequality  is  evident  from  the  fact  that  Y  is  a  feasible 
point  of  P;  the  left  inequality  follows  from  (10). 


t 
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Two  Issues 

Thus  far  we  have  left  two  questions  unanswered;  they  are: 

(i)  How  does  one  cope  with  the  last  segment  of  the  linear  spline  if 
it  generates  an  error  less  than  e,  and  how  does  one  use  such  a  fact 
to  refine  the  approximation?  (ii)  How  would  one  proceed  if  one 
desires  to  weigh  the  various  activities  differently,  as  reflected 
in  the  permissible  error  e^? 

Rather  than  clutter  this  paper  with  a  mass  of  symbolism  and 
formulas  that  will  undoubtedly  mask  the  rather  elementary  nature  of 
the  response  to  these  two  issues,  we  prefer  to  respond  within  the 
context  of  a  numerical  example  that  is  deliberately  chosen  for  its 
simplicity  and  tractability.  The  reader  should  have  no  difficulty 
in  extrapolating  to  more  complex  functions.  Indeed,  many  of  our 
formulas  developed  for  this  simple  example  are  directly  applicable 
to  others,  as  wc  shall  demonstrate. 

NUMERICAL  EXAMPLE 

Perhaps  the  simplest  nonlinear  convex  decreasing  cost  function 

is  a  quadratic  with  continuous  derivative  over  fO,  »]  of  the  form 
2 

a  +  B(u  -  v)  .  Consider  a  project  network  with  such  activity  cost 
and  other  data  given  in  Table  1  and  Fig.  3. 


Table  l  approximately  here 


/ . 

-V  C-- 
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Figure  3  approximately  here 

This  problem  was  solved  optimally  by  the  methods  of  [  ];  the 

optimal  compression  function  is  shown  in  Fig.  4. 

Figure  4  approximately  here 

Now  it  is  desired,  for  purposes  of  illustration,  to  solve  the 

problem  by  approximating  the  cost  function  Cj(y^)  =  (u^  - 

2 

y^.)  using  the  algorithm  for  finding  the  best  FDIS.  The  error  e 
is  to  be  chosen  to  guarantee  that  the  optimum  value  of  the  approxi¬ 
mating  program  does  not  deviate  from  the  optimum  value  of  the  original 
problem  by  more  than  10%  of  the  true  optimum,  i.e.,  we  are  adopting 
Criterion  2  with  w  =  0.10. 

It  is  easy  to  see  that  the  graphical  construction  of  Fig.  2 
can  be  translated  into  the  following  analytical  steps: 

1.  Determine  y  (the  point  of  maximum  error  in  the  kth  segment  of 
the  approximation)  from 

[C (vk)  ♦  e  -  C(yk_1)]/[yk  -  yk_1]  =  C'(yk)  ; 

k  =  1,2,  .  .  .  (14) 

where  C'(yk)  =  dC/dy|^_  -k.  For  the  case  C(y)  =  a  +  B(u  -  y)“, 
the  above  expression  yields 
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yk  =  yk  1  +  /  c/fi  ;  y°  =  «.  (14') 

and  the  slope  of  the  line  segment  between  yk  1  and  yk  is  given 
by 

C'(yk)  =  -  2B(u  -  yk_1  -  /e7¥)  (15) 

2.  Determine  y  (the  end  point  of  the  kth  segment)  from 

[C (yk)  -  C(yk_1))/[yk  -  yk_1]  =  C'(yk)  (16) 

For  the  above  quadratic,  this  yields 

yk  =  yk_1  +  2  /m  (H>’) 

Finally,  using  (12)  with  w  =  0.10,  m  =  5  and  C(U)  =  11  (=^<*.j) 
we  have  e  =  0.220. 

Substituting  for  this  value  of  e  and  the  other  network  parame¬ 
ters,  we  immediately  deduce  the  linear  splines  given  in  Table  2. 


Table  2  approximately  here 


To  help  interpret  this  table,  consider  arc  (2,3):  The  spline 
has  four  segments  as  shown  in  Fig.  5,  with  the  given  knots  and  slopes. 
The  variable  y2  3  is  now  replaced  by  the  four  variables  Xj2>  x13,  Xj4 
and  x15,  whose  coefficients  in  the  (linear)  criterion  are  precisely 
the  corresponding  "slopes"  of  Fig.  5. 
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Figure  5  approximately  here 

Program  P  in  the  original  variables  {y^}  is  now  translated 
into  the  approximating  linear  Program  S: 

Minimize  H(Ts)  =  93  -  7.062  x:  -  5.187  x2  -  3.310  x3 


-  1 

.434 

X4  - 

.25  x$ 

-  • 

9062 

x,  -  7.186  x- 
n  7 

-  5 

.310 

x8 

3.434 

X9 

1.560  x10  -  .323  x 

-  5 

.062 

X12 

-  3.186 

X13 

-  1 

.310  x,. 

14 

—  • 

188  x15  - 

7.062 

X16 

-  5. 

187  xj7 

-  3 

.310 

X18 

-  1.434 

X19 

—  • 

25  X20 

-  7 

.062 

X21 

-  5.187 

x22 

-  3 

.3!0  x23 

-  1 

.434 

x24 

-  .25  x 

25 

Subject  to: 


*1 

'  *2  +  X1  +  X2  +  X3 

+  X4  +  XS 

< 

0 

*3  + 

X6  +  X7  4  X8  +  X9  + 

X10  +  X1 1 

< 

0 

t 

2  ”  t3  +  X12  X13  + 

X14  +  X15 

< 

0 

t2 

+  x16  4  x17  +  x18  4 

X19  4  X20 

< 

T 

s 

*3 

+  X21  +  X22  +  X23  + 

X24  4  X25 

< 

T 

s 

X1  +  X2  +  X3 

4  X4  4  X5 

< 

4 

X6  +  X7  +  x8  +  X9  + 

xio  4  X11 

< 

5 

14 


Xj.  _<  .248; 

x.  <  .938, 
i  — 


+ 

12 

X13  + 

X14  + 

X15 

< 

3 

+ 

x  ^  + 

x,rt  ♦ 

x.A 

< 

4 

17 

18 

19 

20 

— 

^  + 

x  + 

x  + 

x 

< 

4 

22 

23 

24 

25 

x12  _<  .31;  Xj  j-  <_  .186;  xof)  <_  .248;  Xj,.  .248,  and 

for  all  other  i;  x.  >  0,  fcr  all  1. 

’  i  — 


This  LP  was  solved  for  the  three  values  of  Ts  given  by  the 
"breakpoints"  in  the  optimal  compression  function  of  Fig.  4,  namely 

Tg  =  8,  6.75  and  4.  The  results  are  given  in  Table  3,  which  contains 

* 

one  additional  result  for  comparison,  namely  the  column  labeled  L  (Ti_). 


Table  3  approximately  here 


This  column  gives  the  optimal  solution  under  the  "crude"  linear 
approximation  of  each  function  by  the  single  line  segment  joining 
the  two  end  points  of  C.^(y^j)  (see,  e.g..  Fig.  5).  The  values  of 
the  absolute  error  between  the  "true"  and  "approximate"  optima  are 
given  in  column  5;  and  their  relative  magnitudes  (to  the  "true" 
optima)  are  given  in  the  last  column.  The  latter  are  seen  to  be  well 
within  the  0.10  specified  tolerable  error.  Their  small  values  are  a 
reflection  of  the  conservatism  induced  by  having  bared  the  value  of 
e  on  the  lower  bound  C(U);  see  Eq.  (12). 
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Table  3  contains  one  additional  piece  of  information,  namely 

* 

the  a  posteriori  bounds  on  the  optimum  C  (y);  see  liq.  (13).  The 
bounds  are  rather  good,  and  they  do  indeed  contain  the  "true"  value. 

It  is  also  instructive  to  compare  the  optimal  activity  dura¬ 
tions  under  the  three  representations  of  C-j(yjj):  The  original 
quadratic  function,  the  FOIS,  and  the  "crude"  linear  approximation. 
This  is  given  in  Table  4.^  It  is  evident  that  the  differences  among 
the  activity  durations  remain  significant  throughout  the  ranges  of 
T  between  the  "Quadr."  and  "Linear  Spline"  on  one  hand  and  the 
"Crude  Linear"  on  the  other,  with  the  difference  between  the  first 
two  being  rather  small. 


Table  4  approximately  here 


We  are  now  ready  to  respond  to  the  two  issues  raised  above. 

The  Last  Segment  and  Improved  Approximation 

First,  there  is  the  question  of  the  "last  linear  segment"  of 
the  spline  that  terminates  at  u.  A  glance  at  Fig.  5,  for  example, 
and  some  simple  calculations  reveal  that  this  last  segment  (in  the 
interval  [2.814,  3])  has  a  maximum  error  of  only  .008649,  a  far  cry 


+The  fourth  column  under  each  T  ,  labeled  "Improved  Linear  Spline", 
will  be  explained  below. 


* 
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from  the  allowable  value  of  e  =  0.220.  Consequently,  we  may  obtain 
a  smaller  error  c  throughout  the  range  of  y  if  we  insist  on  adding 
one  more  restriction  to  the  FDIS;  namely,  that  it  terminates  on  the 
function  C..(y..)  at  both  l.  .  and  u.  .  while  utilizing  the  same  number 
of  knots.  The  existence  of  such  spline  is  guaranteed  by  the  finite¬ 
ness  of  the  range  of  y.^'s. 

Fortunately,  for  the  quadratic  cost  function  adopted  in  this 
example,  the  analytical  determination  of  the  optimum  FDIS  satisfying 
this  additional  condition  is  rather  elementary.  (This  was  another 
reason  for  the  adoption  of  this  particular  cost  function.) 

Consider  the  first  segment  of  the  FDIS  between  y°  and  y^ .  It 
is  easy  to  see  that  (see  Fig.  3): 


1  =  LCzVlJLcA  re-1  -  v°i  - 


1  0 

v  -  v 


(y 1  -  yv)  -  [c(yx)  -  c(y°)] 


v°  =  1 


(17) 


But,  since  C(y)  -  a  +  8(u  -  y)“,  we  have 


rt  i  i  ,  0.  ,  0  lw_  0  1. 

C(y  )  -  (. (y  )  =  6(y  -  y  )(2u  -  y  -  y  ) 


and  similarly. 


C(y  S  -  C(y°)  =  B(y°  -  y1)  (2u  -  y°  -  y1) 


Substituting  in  the  expression  for  c*,  we  obtain 


c1  =  B(y°  -  y*) (y1  -  y1) 


'1 
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But  y  is  the  point  at  which  the  tangent  to  C(y)  is  parallel  to  the 
line  segment;  i.e.. 


dC(y) 

dy 


-1 

y=y 


=  -  2g(u  - 


y1) 


=  £CQ  - 


c(y°) 


y  -  y 


=  -  B(2u  -  yn  -  y1)  (18) 

yielding  the  relation  between  y*  and  y*  (recall  that  y^  is  fixed 
at  S.) : 

y1  =  (y°  +  yX)/2 


Consecjuently, 


1 

e 


2 


By  a  similar  argument,  we  can  easily  deduce  that,  in  general, 
t  k  k-1 

e  =  g (?— Y - )  1  .  .  .  ,n+l 

Now  we  impose  the  desired  conditions,  which  are  three:  (i)  =  c, 

a  constant  for  all  k;  (ii)  y°  =  l;  and  (iii)  yn+1  =  u.  They  com¬ 
pletely  determine  the  value  of  e.  We  proceed  recursively: 

c1  =  B(^-^)2  =  e  y1  =  2  /?7F  *  l 

e2  =  g(^--~  ^-)2  =  e  y2  =  2  /U&  +  y1 

=  2x2  /e/g  ♦  Jt 


A 
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and  in  general  =  B(^ - ^ - )“  y^  =  k  x  2  /e/B  +  S. 

Finally,  we  have  vn  *  =  u  =  (n+1)  x  2  /?/s  +  l 


which  yields 


e 


8  [■ 


u  -  S., 2 
2(n+l)  ' 


(10) 


This  is  the  new  value  of  e  to  be  used  in  the  construction  of  the 
FDIS  for  this  activity.  Note  that  now  the  maximum  error  of  approxi¬ 
mation  may  be  different  for  different  activities. 

To  illustrate,  consider  once  more  activity  (2,3)  of  Table  1 
and  Fig.  5;  We  have  n  =  3  knots,  6  =  1,  u  =  3  and  8.  =  0;  hence 

9 

e  =  (3/8)  =  9/64  -  0.140625  (compare  with  the  previous  value  of 
e  =  0.220).  With  this  value  of  c,  straightforward  substitution  in 
Fqs.  (14),  (15)  and  (16)  yields  the  following  Table  5: 


Table  5  approximately  here 


The  new  knots  and  slopes  for  the  various  FOIS's  of  the  project 
of  Table  1  are  given  in  Table  (>.  Of  course,  the  number  of  x-variables 
for  each  activity  remains  the  same  as  before,  though  their  ranges  anil 
coefficients  in  the  linear  criterion  are  now  different.  The  new 


Table  (>  approximately  here 
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program,  denoted  by  S,  was  solved  for  the  same  three  values  of  Tg 
as  before,  namely,  Tg  =  8,  6.75  and  4.00,  and  the  optimum  values  are 
summarized  in  Table  7,  together  with  the  corresponding  values  for 
the  "true"  optimum  and  for  the  "original"  FDIS.  As  to  he  expected, 
the  absolute  and  relative  errors  are  smaller  than  before.  Moreover, 
the  bounds  on  the  optimum  (also  given  in  Table  7)  are  much  sharper. 
(For  comparison,  see  Table  3.)  The  optimum  durations  of  the  activi¬ 
ties  at  these  various  target  completion  times  have  been  included  in 
Table  4  to  facilitate  comparison  with  other  conditions. 


Table  7  approximately  here 


The  development  thus  far  responds  completely  to  the  question 
raised  concerning  the  "last  segment  of  the  FDIS"  when  the  individual 
activity  cost  function  is  convex  decreasing,  quadratic  and  with 
continuous  derivative  for  yc[£,°°].  Extension  of  the  approach  to 
quadratic  cost  functions  that  possess  discontinuous  derivative  at 
y  =  u  (the  case  treated  in  ref.  [2  ])  is  straightforward  and  will  not 
be  elaborated  upon  here. 

Such  is  not  the  case,  however,  for  general  cost  functions 
that  are  convex  decreasing  in  the  interval  [i,  u],  such  as  the 
function  C(y)  of  (1).  The  reason  is  that  now  we  may  have  to  deal 
with  nonlinear  equations  to  determine  the  new  value  of  c,  though  the 
logic  of  derivation  implied  by  Eqs.  (17)  and  (18)  remains  applicable 
in  toto.  In  such  instances,  the  application  of  simple  iterative 
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schemes  yields  the  desired  result  in  a  few  iterations.  Indeed, 
bounds  on  the  locations  of  the  new  knots  are  easily  deter¬ 

mined  by  applying  the  graphical  (or  analytical)  construction  of 
Eqs.  (14)-(16)  starting  from  the  point  (u,C(u))  rather  than  the 
point  (Jl,C («.)). 

The  concepts  of  iteration  are  best  illustrated  by  example. 
Consider  again  the  simple  network  of  Fig.  3,  and  suppose  that  the 
activity  cost  functions  are  given  by  Eq.  (1).  Elementary  calcula¬ 
tions  (see  also  ref.  [6  ],  pp.  219-22S)  following  the  logic  of 
Eqs.  ( 14) - (16)  yield: 

C(ykl-.;  C  -  ctyb-1)  .  c, ,-k,  ,  „k 
yk  -  yK 

Substituting  for  C(y)  =  a/(b  +  ry)  in  the  above  equation  yields 
a/_(b_^ryk)  +  =  .  ar/(b  +  r^2 

y  -  y  ’ 

or 

[e  -  a/(b  +  ryk_1)]  (b  +  ryk)~  +  a/(b  +  rvk) 

,-k  k-1 

=  -ar  (y  -  y  ) 

-k 

This  is  a  quadratic  in  (b  +  ry  )  whose  solution  yields 

1  /° 

-k  {  f  2  ..  k-1.  ,  , ,,  k-K,  % 

y  =  [-  a  -  (a  +  a (b  +  ry  )  (e  -  a/(h  +  ry  )}  ]  1 

r(e  -  a/(b  +  ry^"*))  -  b/r 
-k  k-1 

(choosing  the  root  that  yields  y  >  y  ). 


/ 
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Determine 

mk  =  C' (y) |y_yk  =  -  ar/(b  +  ryk)2 

and  compute  yk  from 

k  ,  k,  k-K,,  k  k-1. 

y  =  -  (ar  +  m  b  (b  +  ry  ))/m  r  (b  +  ry  ) 

k  k  “  l 

If  y  u  stop  and  the  last  segment  of  the  FDIS  joins  C(y  )  5 

C(yn)  and  C(u),  i.e.,  put  yn+1  =  u. 

To  illustrate,  assume  the  parameters  are  as  given  in  Table  8, 
which  is  constructed  so  that  the  cost  function  satisfies  the  (rea¬ 
sonable)  condition  that  would  permit  comparison  later  on  with 

previous  results,  namely,  that  C(y)  passes  by  the  points  C(£)  and 

2 

C(u)  of  the  quadratic  cost  function,  a  *  g(u  -  y)  .  It  is  easy  to 
see  that  the  cost  function  C(y)  =  a/ (b  +  ry)  may  be  written,  for 
simplicity,  as  C(y)  =  a'/(l  +  r'y)  when  both  numerator  and  denominator 
are  divided  by  b.  Since  =  0  and  6=1  for  all  activities  in  this 
example,  this  condition  yields: 

a'  =  C (0)  and  r'  =  —  -  I],  resulting  in 

u  C(u) 

a  =  aC(0)  ;  b  =  a  ;  r  =  u  (19) 


Table  8  approximately  here 
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For  notational  simplicity,  we  refer  to  the  knots  and  errors 
that  result  from  starting  the  spline  at  the  point  (l,C  (J.) )  as  result¬ 
ing  from  forward  iteration,  and  to  the  knots  and  errors  that  result 
from  starting  the  spline  at  the  point  (u,C(u))  as  resulting  from 
backward  iteration. 

For  illustration,  consider  once  more  activity  (2,3): 

C  (y)  =  1  +  (3  -  y  )2,  which  yields  (see  (19)  above)  C(v)  = 

'<  C.  y  O 

10/(1  +  3y) .  Forward  iteration  is  started  with  the  original  error 
of  Cp  =  0.22,  which  results  in  five  segments  as  shown  in  Table  9. 
(Compare  with  only  four  segments  of  the  quadratic  cost  function.) 


Table  9  approximately  here 


Two  remarks  are  relevant  to  Table  9.  First,  the  unrestricted 
last  segment  of  the  FDIS  "overshoots"  the  upper  bound  on  the  activity 
duration  by  a  considerable  amount.  Therefore,  it  must  be  truncated, 
resulting  in  a  much  smaller  error  than  c^. 

Second,  it  is  rather  easy  at  this  juncture  to  derive  bounds 
on  the  locations  of  the  knots  through  backward  iteration  [i.e.,  start¬ 
ing  with  spline  at  the  point  (u,C(u))  and  proceeding  "backwards"  to 
the  point  (£,C(£))].  This  is  accomplished  in  Table  10  for  the  same 


value  of  eQ  =  0.22. 
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Table  10  approximately  here 


The  bounds  on  the  location  of  the  improved  knots  are  shown  in  Fig.  6. 


Figure  6  approximately  here 


The  easiest  iterative  scheme  is  to  put  the  new  value  of  the  error 
equal  to  the  average  error  in  the  previous  iteration.^  Therefore, 

v  k 

put  e  =  e  /5  =  0.192372.  The  result  is  an  improved  spline 
k 

(still  of  five  segments)  as  shown  in  Table  11. 


Table  11  approximately  here 


V  k 

We  iterate  once  more  with  s  =  /  ,  e  /5  =  0.187169,  and  obtain  the 

2  k  1 


values  shown  in  Table  12. 


^Alternative  schemes,  such  as  choosing  the  end-points  of  the 
just-determined  intervals  of  the  knots,  were  not  systematically 
investigated. 
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Table  12  approximately  here 


The  deviation  from  in  the  last  segment  is  0.000755  or  approxi¬ 
mately  0.4%,  and  iteration  may  very  well  stop  here.  (In  a  computer 
program,  one  would  pre-specifv  an  "acceptable"  deviation  in  the  error 
at  which  iteration  is  terminated.)  We  decided  to  do  one  more  itera- 

tion  with  c,  =  =  0.  187018.  The  results  are  shown  in  Table  13. 

'  ‘k  “ 


Table  13  approximately  here 


To  be  sure,  the  final  knots  are  within  their  respective  bounds  estab¬ 
lished  by  the  forward  and  backward  iterations;  see  Fig.  6.  (For 
example,  y1  =  . 1 13978e[0. 069616,  0.126214),  y2  =  0.298232Gf0. 241950, 
0.340470],  etc.) 

The  complete  (improved)  FDIS's  for  all  five  activities  are 

given  in  Table  14,  obtained  by  similar  (forward)  iterations  over  the 

values  of  the  error  e...  To  be  sure,  the  final  errors  achieved  are 

il 

all  smaller  than  the  starting  error  of  =  0.22,  and  they  differ 
from  activity  to  activity.  The  optimal  solutions  at  the  three 
"breakpoints"  are  given  in  Table  15,  together  with  the  respective 
bounds  on  the  "true"  optimum  values  based  on  them  and  on  the  value 

V 

/  .  e..  =  0.920068. 

(ij)  1J 
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Table  14  approximately  here 


Table  15  approximately  here 


The  following  remarks  are  pertinent. 

It  is  evident  that  the  cost  of  project  compression  is  quite 

sensitive  to  the  form  of  the  assumed  activity  time-cost  trade-off 

function,  contrary  to  the  common  lore.  For  instance  ,  at  T  =  6.75, 

the  cost  of  the  project  varies  from  a  low  of  13.84  for  the  FDIS  to 

the  (steeply)  convex  cost  function  of  (1);  to  17.44  for  the  quadratic 

2 

convex  cost  function  a  +  6(u  -  y)  ;  to  27.70  for  the  FDIS  to  it;  to 
17.85  for  its  "improved"  FDIS;  to,  finally,  27.25  for  the  "crude" 
linear  approximation! 

Though  we  have  not  conducted  the  requisite  empirical  investi¬ 
gation,  it  is  our  belief,  on  the  basis  of  limited  evidence,  that  the 
ratio  of  cost  between  the  "crude"  linear  approximation  and  that  of 
Eq.  (1)  will  be  significant  (  ^  1.25).  (Of  course,  this  depends  on 
the  amount  of  compression  required.  In  the  above  example,  the  dura¬ 
tion  was  reduced  by  approximately  40%. ) 

The  Differential  Weighting  of  Activities 

We  remind  the  reader  that  this  is  the  second  question  raised 
under  "Issues"  above.  It  is  possible,  in  some  application,  that 


«•' 

A  . 
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large  errors  will  be  tolerated  in  some  activities  but  not  in  others. 
This  differentiation  among  activities  may  be  implemented  by  specify¬ 
ing  a  "tolerable"  ratio  between  any  pair  of  errors.  This  will 
introduce  a  slight  modification  to  the  development  of  hqs.  (10) 

through  (12).  In  particular,  we  now  write  e. .  =  a.  .e,  a. .  >  1  and 

— 

let  A  =  /  .  a. ..  Then  we  have 

C - >  -1  i 


(ij) 


ij 


v(S)/v(P)  <  1  +  Ae/v(p)  ;  v(P)  >  0 


which  results  in 


e  oj.C(U)/A  and  =  a^ ^  e  (20) 

where,  as  before,  C(tl)  =  /  .  C..( u..l. 

fit)  1-1 

As  illustration,  consider  the  simple  network  of  Fig.  3  which 
has  served  as  our  vehicle  thus  far,  and  suppose  that  the  weighting 
of  the  activities  is  in  the  ratio:  Yj  2  :  yl  3  :  y2  3  !  y2  4  :  y3  4 

=  1  :  2  :  1  :  1  :  3.  Then  A  =  8  and,  for  the  parameters  given  in 
Table  1  and  with  w  =  0.10,  we  obtain  e  =  (0.1  x  11)  -4-  8  =  0.1375. 
Whence  ^  =  0.1375  =  c?  ^  ^  =  0.2650  and  ^  =  0.4125. 

From  this  point  onwards,  the  procedure  of  spline  approximation  follows 
identical  lines  to  those  of  either  the  FDIS  or  the  improved  FDIS, 
and  will  not  be  repeated. 
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Figure  4  Optimal  Compression  Function 
Under  Quadratic  Cost 
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Table  5  Improved  FDIS  for  Activity  (2,3) 
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1 

0 

C.  75 

0.375 

-5.25 

2 

0.75 

1.50 

1.125 

-3.750 

3 

1.50 

2.25 

1  .875 

-2.250 

4 

2.25 

3.00 

2.625 

-0.750 

(=  U2 , 3 

) 

M 

A 


fc: 


17.0525  <  C  (6.75)  <  17.5050 


Table  8  Parameters  of  Network  Under 
C(y)  =  a/(b  +  ry) 


Activity 


i: 


u .  . 

ID 


a .  . 

ID 


b.  . 
ID 


(1.2)  0  4 

(1.3)  0  5 

(2.3)  0  3 

(2.4)  0  4 

(3.4)  0  4 


57  3 

116  4 

10  1 

36  2 

17  1 


Table  9  Forward  Spline  for  =  0.22, 
Activity  (2,3) 


k 

JS. 

y 

— JS. 

y 

m 

e 

0 

0.00 

- 

- 

- 

1 

0.126214 

0.058052 

-21.760562 

c0 

2 

0.340470 

0.223124 

-10.765040 

c0 

3 

0.748720 

0.520535 

-  4.571900 

E0 

4 

1.681883 

1.143343 

-  1.528651 

£o 

5 

3.00 

(4.659654) 

2.258459 

-  0.496225 

0.0816 

y  :  kth  knot  value 

__k 

y  :  kth  value  of  point  of  max.  error 

Ir  _ Ir 

m  :  slope  C'(y)  evaluated  at  y 

(The  last  number  in  parenthesis  is  the  original  knot  if 
the  value  of  y  were  not  bounded.) 


* 


Table  10 

Backward  Spline 
Activity  (2,3) 

with  Cq  *  0.22 

k 

-k 

Y 

y 

3.00 

- 

- 

1.211247 

1.935700 

-0.647425 

0.554096 

0.837438 

-2.431836 

0.241950 

0.381176 

-6.529248 

0.069616 

0.148133 

-14 . 379588 

0.00 

0.013118 

-27.771252 

(-0.035459) 


Table  11  Forward  Spline  for  c ^  =  0.192373, 
Activity  (2,3) 


k 

y 

1  >, 

k 

m 

k 

e 

0.00 

- 

- 

- 

0.115999 

0.053678 

-22.255222 

£i 

0.305046 

0.202245 

-11.620695 

ei 

0.644340 

0.456683 

-  5.340799 

ei 

1. 348390 

0.948921 

-  2.027357 

E  i 

3.00 

2.034312 

-  0.594628 

0.166357 

(3.214741) 

A 


Table  13 

Forward  Spline  for 
Activity  (2,3) 

t3  =  0.187018, 

k 

-k 

k 

k 

k 

,  y 

i 

y 

m 

c 

0 

l 

0.00 

- 

- 

- 

1 

0.113978 

0.052806 

-22.355783 

e  3 

2 

0 . 298232 

0.198180 

-11.799138 

c  3 

3 

0.625102 

| 

0.444686 

-  5.506778 

E  3 

4 

1.291157 

0.914452 

-  2.140912 

e  3 

5  1 

3.00 

1.993678 

-  0.615576 

E  o 

(3.000008) 


(Deviation  <  10 


Table  15  FDTS  Approximation  for  C(y)  =  a/(b+ry) 


T  =  8  1 

C 

T  =6.75 
s 

o 

o 

It 

-  W 

FDIS  Cost 

12.5729 

13.8444 

20 . 7589 

y  1  <  2 

‘  4.0 

2 . 750 

1 . 8600 

Yl,  3 

1 

!  5.0  i 

i 

4 . 96784 

3.1512 

Y2,  3 

2.21784  1 

i  i 

2.21874 

I 

|  1.2912 

V2 , 4 

1  4.0  i 

,  1 

1  ! 

4.00 

/ 

|  2.1 39^ 

i 

Y  3 , 4 

1 

1  1.78216  I 

!  i 

1.78216 

i 

i  0.84876 

11.6528  < 

C*  (  8)  < 

12.39  80 

12.9243  <  C*  (  6 . 7 5)  <_  13.4917 

<■  C* ( 4 . 0 )  <  21.6^42 


19.8388 
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